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ABSTRACT 

Particle theorists typically use expectation values to study the quantum back- 
reaction on inflation, whereas many cosmologists stress the stochastic nature 
of the process. While expectation values certainly give misleading results for 
some things, such as the stress tensor, we argue that operators exist for which 
there is no essential problem. We quantify this by examining the stochastic 
properties of a noninteracting, massless, minimally coupled scalar on a locally 
de Sitter background. The square of the stochastic realization of this fleld 
seems to provide an example of great relevance for which expectation values 
are not misleading. We also examine the frequently expressed concern that 
signiflcant back-reaction from expectation values necessarily implies large 
stochastic fluctuations between nearby spatial points. Rather than viewing 
the stochastic formalism in opposition to expectation values, we argue that it 
provides a marvelously simple way of capturing the leading infrared logarithm 
corrections to the latter, as advocated by Starobinsky. 
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1 Introduction 



Schrodinger was the first to suggest that spacetime expansion can lead to 
particle production by ripping virtual particles out of the vacuum [1]. Fol- 
lowing early work by Imamura |2] , the first quantitative results were obtained 
by Parker [3]. He found that the effect is maximized during accelerated ex- 
pansion, and for massless particles which are not conformally invariant [1], 
such as massless, minimally coupled scalars and (as noted by Grishchuk [5]) 
gravitons. Precisely this process is responsible for the primordial spectra of 
scalar and tensor perturbations which are believed to arise from inflation [6] , 
the scalar contribution of which has been imaged [7]. 

Inflationary particle production results from the background gravitational 
field acting on quantum matter (and graviton) fluctuations. It is natural 
to wonder about the complementary process of back-reaction in which the 
newly produced particles modify the background gravitational field, either 
directly or through their self-interactions. People who approach this from the 
perspective of particle physics typically attempt to quantify back-reaction 
using expectation values or in-out matrix elements [H [9l [TOl [TT|, |12l [131 El 
[T5l[l6lI7lE8l[T9]. 

Many cosmologists dismiss the use of expectation values and in-out ma- 
trix elements as giving an unreliable average over vastly different portions 
of a quantum wave function which has actually decoheredQ They fear that 
particle theorists are falling victim to a sort of cosmological Schrodinger Cat 
Paradox based on a fictitious, mean geometry which bears no relation to 
what any observer would experience. Cosmologists prefer to instead study 
back-reaction using a stochastic formalism in which the super-horizon modes 
of various fields are regarded as classical, random variables pTl [2^ [2 H [25] . 
Whereas expectation values in a homogeneous and isotropic state necessar- 
ily produce a homogeneous and isotropic geometry, cosmologists assert that 
the actual universe is not even approximately homogeneous on super-horizon 
scales [23]. They also fear that if the back-reaction inferred from expecta- 
tion values ever becomes significant then the resulting universe would show 
unacceptable spatial fluctuation. 

There is no question that cosmologists are right about certain operators 
being poorly described by their expectation values. For example, the vacuum 

^For an excellent recent study of cosmological decoherence, which does not necessarily 
endorse the anti-VEV position, see |20j . 
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expectation value of the stress-energy tensor is homogeneous and isotropic, 
whereas we perceive inhomogeneities and anisotropics. What isn't clear is 
whether or not expectation values are in all cases misleading. Such an ex- 
treme position would be embarrassing for cosmologists (in spite of the fact 
that some do advocate it) because the primordial power spectra are defined 
by taking expectation values [25] . 

We suspect that the reliability of expectation values depends upon the 
operator under study. For operators which average to zero, such as the den- 
sity perturbation, the entire result arises from the decoherence effect, so one 
makes an enormous mistake by ignoring it. Other operators — for example, 
the square of a scalar field — acquire a significant homogeneous expectation 
value upon which spatial variations are superimposed. Any quantum fluctu- 
ation drives this sort of operator positive, so one might happen to inhabit a 
special region of the universe in which there is little effect for a long time, 
but there will sooner or later be a large effect. The expectation value of such 
an operator can correctly reflect the long-term trend everywhere in space, 
even though it misses variations from one region to another. 

The purpose of this paper is to use the noninteracting, massless, mini- 
mally coupled scalar on de Sitter background to give a quantitative assess- 
ment of the two key issues under dispute between particle theorists and 
cosmologists: 

1. How unreliable are expectation values? and 

2. How much spatial variation should one expect? 

In section 2 we review the de Sitter geometry, the scalar field model, and 
our stochastic realization of it. In section 3 we demonstrate that the square 
of the stochastic field follows a distribution whose mean grows with the 
number of e-foldings. Although fluctuations about this mean are significant, 
they do not contradict the picture provided by expectation values. In section 
4 we address the issue of spatial fluctuations by studying the difference of 
the stochastic scalar at two points held at a fixed physical distance from one 
another. Our conclusions comprise section 5. 

2 The Theoretical Context 

We work on the spacetime manifold x R. We assume each of the three 
spatial coordinates lies in the range, < < where H is the Hubble 
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constant. We also assume that the geometry is locally de Sitter, 



ds^ = —dt^ + a^{t)dx ■ dx 



a{t) 



(1) 



Note that the scale factor a{t) is normalized to unity at t = 0, rather than 
at the current time. Because the geometry ([T]) is homogeneous and isotropic 
we can expand any function of spacetime in a spatial Fourier series, 



ik-x 



k = 2TiHn 



(2) 



A key event for any mode of wave number k = \\k\\ is the time tk of hori- 
zon crossing when its physical wave number redshifts down to the Hubble 
constant, 

-^ = H. (3) 



The Lagrangian for a massless, minimally coupled scalar is, 
C 



(4) 



For all nonzero wave numbers k the canonically normalized, Bunch-Davies 
mode functions are. 



u{t, k) 



H 



1 



ik 



Ha{t) 



exp 



ik 



Ha{t) 



uu — uu 



(5) 



For the = mode the two solutions are a constant and l/a^{t). Hence the 
field operator can be expanded as. 



Lp{t,x) = H^\Q 



u{t, k)e"^-'^AH + u*{t, k)e-''^-^Al 



where the nonvanishing commutators are 



A ^ 4 ' 



(6) 



(7) 



The state \ which is annihilated by all is known as Bunch-Davies vac- 
uum. It does not matter very much what we assume about its dependence 
upon the 0-mode coordinate. 
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The quantity f{t,x) is a quantum field operator which obeys the Uncer- 
tainty Principle, 

_ip{t,x),ip{t,y)]=i6'^{x-y) . (8) 

Because the free field expansion (|6]) contains arbitrarily large wave numbers, 
expectation values of coincident products of '^(t,x) can harbor ultraviolet 
divergences. All of these features are absent in the stochastic realization of 
ip{t, x) which we construct by taking the infrared limit of the mode functions, 

lim uit, k) = , (9) 

and retaining only the super-horizon modes. We denote this stochastic field 
as (p{t,x) and its definition is [17], 



Instead of being creation and annihilation operators, and At are consid- 
ered to be (complex conjugate) stochastic random variables which follow a 
normal distribution with mean zero and standard deviation one. That is, 
one can express as the sum of two real, independent random variables, 

Afi = afi + iPri where p(a^ = x,Pfi = y) = . (n) 

3 Fluctuations about the Mean 

It is useful to express if{t,x) in terms of the real stochastic variables and 
/3ri which were introduced in ( ITTj) . 



(p{t, x) = ^/ -^j-^ 9{Ha{t) — k] cos{k-x) afi — sm{k-x) . (12) 



n^O 



In this form one can recognize (p{t, x) as the sum of a vast number N{t) of 
independent Gaussian random variables, 

/ \ /•a/27r a'^(t) 

N{t) = 2 ^ d[Ha{t)-k) ~ 2 X 47r dnn^ ^ . (13) 

-./rv 1 /27r StT 



Of course the sum of any number of independent Gaussian random variables 

gives another Gaussian random variable whose mean is the sum of the means 
and whose variance is the sum of the variances. Because the mean of each 
variable is zero, the mean of (p{t, x) vanishes. On the other hand, its variance 
grows, 

= Y.7TU^{Ha(t)-k)\cos\k-x)+sm''{k-x)] , (14) 

pa/2TT }J5 J^2 

~ 47r / dnn^ — — — — = — \n[ait)] . (15) 

Hence we can say that, at any given spacetime point, the stochastic field 
(p{t,x) follows a normal distribution with mean zero and variance cr^(t), 

''(^"•""' = '') = 7iw"'''["^)- '''' 

The fields at different spacetime points are not statistically independent, but 
we will not need to worry about that until the next section. 

The variable 0{t,x) provides a classic example of why cosmologists dis- 
trust expectation values. Although its expectation value vanishes, the sto- 
chastic field experiences very significant and growing fluctuations, as its vari- 
ance reveals. Someone interested in the behavior of (p{t,x) would derive a 
completely misleading picture from its expectation value. 

On the other hand, consider the variable (p'^{t, x). Because it is the square 
of a Gaussian with mean zero it follows a distribution whose mean is the 
variance of the original Gaussian, 

1 r Z 1 

p((p^(t,x) ^ Z) = , exp — . (17) 

^ ^2na^{t)Z ^ 2a2(t)J ^ ^ 

The variance of (p'^{t,x) might seem to vindicate the extreme cosmologist 
position that expectation values are never reliable, 

(^{0\t,x) - {0\t,x))y^ = //^^= M-^] = • (18) 

However, this only implies that that stochastic fluctuations about the mean 
are significant, even when considered as a fraction of the mean. 
What the expectation value {0^{t,x)) — cr'^{t) really tells us is: 
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• That (p'^(t,x) grows without bound; and 

• That this growth is proportional to the number of e-foldings, ln[a(t)]. 

We can gain a quantitative assessment of the rehabihty of the first conclusion 
by using (fT7|l to compute the probability that (p'^it, x) remains less that some 
constant value 



Prob((^2(t, 



x) < 




(19) 



(20) 



Because cr'^it) grows like the number of e-foldings, we see that the probability 
for 0^{t,x) to fall below any fixed value $^ goes to zero at late times. Of 
course that vindicates the inference of growth without bound. The second 
inference can be tested by computing the probability for (p'^{t,x) to be above 
some time dependent value ^^(t) which grows faster than cr^(t), 



Prob(^2(t,f) > <^^{t) 



oo 1 

dZ 



exp 



$2(t) 



O 



<l>2(t) 



(21) 



(22) 



Under the assumption that o"^(t) / $^(t) goes to zero we see that this probabil- 
ity also approaches zero. Hence the second inference is equally valid, and we 
conclude that no serious error arises from using expectation values to study 

One might object that there is still a substantial disagreement between 
quantum field theoretic expectation values and stochastic samples because 
the former contain an ultraviolet divergent constant in addition to the in- 
frared logarithm [27] . 



n 



(f {t,x) 



n 



Divergent Constant ) + 



47]-2 



ln[a(t)] 



(23) 



The form of this divergence depends upon the regularization technique; with 
dimensional regularization in D spacetime dimensions one finds [28] , 



Y{D-1 

(47r)T 




21n[a(t)]-^(l-^) 
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+^(^) + ij{D-l) + ^(1) + 0(a-2) |. 



(24) 



However, many fully dimensionally regulated and renormalized computations 
have been done involving massless, minimally coupled scalars [291 EOl [311 [321 
[331 [M] and gravitons [351 [SHI [37] with various interactions. What always 
happens is that counterterms absorb the ultraviolet divergence and leave the 
infrared logarithm as the dominant contribution to the final result. That is 
just what one gets from {(p'^(t,x)). 

One might also object that the example of (p'^{t,x) is contrived because 
it represents two fields at the same spacetime point. However, it is exactly 
these terms which are most responsible for the growth of the vacuum energy 
in \(p^ theory [22] and for the photon developing a mass in scalar quantum 
electrodynamics [3B1 [3D] • Note also that the use of expectation values gives 
precisely the correct results for the leading logarithm effects in each case 
[291 [311 [30], 



X 



El 

47r2 



ln[a(t)] g^,, (25) 



Mi 



(^^(t,f) 



47r2 



ln[a(t)] 



(26) 



This is not an accident, nor is the coincidence restricted to lowest order 
perturbative results such as those given above. For scalar potential models 
one can show that Starobinsky's formalism [21j captures the leading infrared 
logarithms of quantum field theoretic expectation values to all orders [39] . 
The best way of viewing the stochastic formalism is not as an alternative 
to expectation values but rather as marvelously simple way of deriving the 
most important contributions to them. 

In some cases the stochastic formalism can do even more. Starobinsky 
and Yokoyama have shown how it can be used to sum the series of leading 
infrared logarithms to derive nonperturbative results for the late time limit 
[iO] . For example, these results explicitly disprove the two simplest imple- 
mentations of the common notion that cosmological evolution can be viewed 
as a renormalization group flow [H]. The Starobinsky- Yokoyama technique 
has recently been extended to Yukawa theory [33] and to scalar quantum 
electrodynamics [12]. It has not yet been extended to gravity but there are 
reasons for believing that some version of it can be [l3l 
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4 Fluctuations in Space 



We can study spatial variation by taking the difference of the fields at points 
on the surface of simultaneity, 



Aif{t, Ax) = (p{t, 0) - (p{t, Ax) 



I'^H^n/rr ,\ • /k-Ax^\ . fk'Ax^ f^'Ax^^ 
—^e\Ha[t)-kj sin(^ ^ j sin(^ ^ j an + cos(^ ^ j /^^ 



(27) 
.(28) 



Just like i^(t,x), this is a sum of independent Gaussians, so Ai^(t,Ax) is 
itself a Gaussian. Because the mean of each constituent is zero, the mean 
of A(p{t, Ax) also vanishes. Its variance is the sum of the variance of each 
constituent. 



r 

1 



sin(/cAx) 



ya/27r 

Ji/2n {2'KHn)^ L kAx 

— - / 

if2 r sm[a{t)HAx] sin[i7Ax] 



(29) 

(30) 
(31) 



27r2 [ a{t)HAx 



HAx 



-ln[a(t)] - ci[a{t)H Ax\^ci[H Ax] \. (32) 



The symbol ci(x) stands for the cosine integral whose definition and expan- 
sion for small x are. 



ci{x) = — dt 



cos(t) 



ln(x) + 7 + X! 



n=l 



ln(x) + 7 + 

-1)V" 
2n-2n! ' 



cos(t) — 1 



(33) 
(34) 



Note the appearance of Euler's constant, 7 ~ 0.577215665. 

Quantum gravitational perturbation theory breaks down when the expec- 
tation value of G(p'^{t,x) grows to be of order one [44j, 



Gip\t,x: 



An"' 



X ln[a(t)] 



ln[a(t)] 



GH^ 



> 1 . (35) 



One might expect that this is also when back-reaction becomes significant. At 
this point the fluctuation of each of the two fields in A(p{t, Ax) is enormous. 
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as one can see from expression fllSp . Cosmologists frequently express the 
worry that, if back-reaction of this sort ever becomes significant, it must 
induce similarly large fluctuations in A(p{t, Ax). That would lead to an 
unacceptable level of inhomogeneity in the post- inflationary universe. 

We can test the cosmologists' fear by choosing Ax so as to keep the 
physical distance a constant fraction K of the Hubble length, 

a(t)Ax = I . (36) 

From expression (!32|) and the asymptotic expansion we see that the 
variance rapidly approaches a not especially large constant, 

a'it,Ka-') = ^ — ^ ^-3^+ln(a)-ci(ir) + ci(ira-i) , (37) 



J sm{K) 
2^1 K 



1 + \n{K) + 7 - ci{K) + 0(a-2) . (38) 



This should be compared with the variance of cr^(t) = H'^/Att'^ x ln[a(t)] of 
each field in the difference. So it is false that spatial fluctuations between 
nearby points become enormous whenever back-reaction is significant. 

Although our result (1381) might seem surprising, there is a very simple 
and general reason for it based upon causality. The stochastic field acquires 
its fiuctuations one instant at a time, as each new complement of modes 
experiences horizon crossing and contributes to the stochastic jitter. The 
result at any spacetime point (t, x) depends upon what happened in the past 
light-cone of that point. After a long period of infiation, two fields a fixed 
distance apart very largely share the same past light-cone, so they experience 
almost the same fiuctuations. This remains true even if the faction K is much 
greater than one, because significant back-reaction requires the staggering 
number of ~ 10® e-foldings. 

Let us also note that the level fl38l) of inhomogeneity we see is about 
right to explain the primordial power spectra [H] . So far from the stochastic 
formalism invalidating theories of infiation which are based on back-reaction, 
it provides an essential ingredient. 



5 Epilogue 

We have employed a very simple scalar model on de Sitter background to ex- 
amine the chief criticisms against using expectation values and in-out matrix 
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elements to infer back-reaction: 

1. That stochastic fluctuations about the mean value change the entire 
picture; and 

2. That significant back-reaction necessarily implies an unacceptable level 
of inhomogeneity. 

Neither criticism is supported by our study. In section 3 we found that 
although there is significant stochastic fluctuation in the quantity <^^(t, x), 
there is zero probability either for this quantity to remain bounded, or for it 
to grow at a faster rate than that predicted by its mean value. In section 4 
we showed that there is no growth in the variance of the difference between 
two stochastic fields held at a fixed physical distance from one another. 

The cosmologists' objection that expectation values give a misleading 
picture for some operators is certainly correct for (p{t,x). On the other 
hand, expectation values give a fair representation of (p'^it, x). So the picture 
that emerges is more complex than either a total refutation of expectation 
values or their complete vindication. The fair conclusion would seem to be 
that the sensitivity of a particular mechanism of back-reaction to stochastic 
fluctuations should always be checked, but there are no grounds for rejecting 
the mechanism prior to such a check. 

We do not view the stochastic formalism of Starobinsky as an alter- 
native to quantum field theory but rather as a wonderful tool for isolating 
the most significant corrections |39]. The really intriguing thing is that the 
stochastic formalism can sometimes be used to obtain nonperturbative re- 
sults |ini ESI 112] • We believe this can be done for quantum gravity [121 SI] 
and that it is worthwhile attempting to anticipate the result [45] . 
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